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Abstract 

We investigate the continuum limit of the low lying static-light meson spectrum using Wilson 
twisted mass lattice QCD with Nf = 2 dynamical quark flavours. We consider three values 
of the lattice spacing a ~ 0.051 fm, 0.064 fm, 0.080 fm and various values of the pion mass in 
the range 280MeV<mps £640MeV. We present results in the continuum limit for light cloud 
angular momentum j = 1/2, 3/2, 5/2 and for parity V = + , —. We extrapolate our results 
to physical quark masses, make predictions regarding the spectrum of B and Bg mesons and 
compare with available experimental results. 



1 Introduction 



A systematic way to study B and Bg mesons from first principles is with lattice QCD. Since 
arrih > 1 at currently available lattice spacings for large volume simulations, one needs to use for 
the b quark a formalism such as Heavy Quark Effective Theory (HQET) [H [2] or Non-Relativistic 
QCD j3]. An alternative procedure has recently been proposed [3] which is based on HQET but 
does not make use of the static point. Here we follow the standard HQET route, which enables 
all sources of systematic error to be controlled. 

In the static limit a heavy-light meson will be the "hydrogen atom" of QCD. Since in this 
limit there are no interactions involving the heavy quark spin, states are doubly degenerate, i.e. 
there is no hyperfine splitting. Therefore, it is common to label static-light mesons by parity 
V and the total angular momentum of the light degrees of freedom j with j = |/ it 1/2|, where 
I and ±1/2 denote respectively angular momentum and spin. An equivalent notation is given 
by l±, which reads S = (1/2)", P_ = (1/2)+, P+ = (3/2)+, D_ = (3/2)", D+ = (5/2)", 
F_ = (5/2)+, = (7/2)+ ... The total angular momentum of a static-light meson is either 
J = j + 1/2 or J = J — 1/2, where both states are of the same mass. Note that in contrast to 
parity, charge conjugation is not a good quantum number, since static-light mesons are made 
from non-identical quarks. 

The static-light meson spectrum has been studied comprehensively by lattice methods in the 
quenched approximation with a rather coarse lattice spacing [5]. Lattice studies with Nj = 2 
flavours of dynamical sea quarks have also explored this spectrum [6l [3 El [U [lOl [HI [12] . Here 
following our initial study [131 [H] , we use Nf = 2 flavours and are able to reach lighter dynamical 
quark masses, which are closer to the physical u/d quark mass, so enabling a more reliable 
extrapolation. Note that in our formalism, maximally twisted mass lattice QCD, mass differences 
in the static-light spectrum are 0{a) improved so that the continuum limit is more readily 
accessible. We now extend our study to include three different lattice spacings, which gives us 
confidence that we are indeed extracting the continuum limit. 

In this paper, we approach the B meson spectrum by concentrating on the unitary sector, where 
valence quarks and sea quarks are of the same mass. This is appropriate for static-light mesons 
with a light quark, which is u or d. 

We also estimate masses of Bg mesons with s quarks of physical mass, where the s quark is 
treated as a valence quark in the sea of light u and d quarks (so this is a partially quenched 
study) . We took our s quark mass values from ETMC studies of strange mesons |15t [16] . 

Within the twisted mass formalism, it is feasible to use Nf = 2 + 1 + 1 flavours of dynamical 
sea quarks, which would give a more appropriate focus on the static-strange meson spectrum if 
strange quark sea effects were significant. This is under study by ETMC. 

In HQET the leading order is just the static limit. The next correction will be of order l/niQ, 
where mg is the mass of the heavy quark. This correction is expected be relatively small for b 
quarks, but larger for c quarks. Lattice methods to evaluate these I/ttiq contributions to the 
B meson hyperfine splittings have been established and tested in quenched studies [171 [TSj [TU 
[20} [21] . We intend to explore these contributions using lattice techniques subsequently. An 
alternative way to predict the spectrum for B and Bs mesons is to interpolate between D and 
Dg states, where the experimental spectrum is rather well known, and the static limit obtained 
by lattice QCD assuming a dependence as l/mg. Thus the splittings among B and Bg mesons 
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should be approximately mc/mb ~ 1/3 of those among the corresponding D and Dg mesons. 

For excited Dg mesons, experiment has shown that some of the states have very narrow decay 
widths [22]. This comes about, since the hadronic transitions to DK and DgM (where M is 
a flavour singlet mesonic system, e.g. r/', vrvr or /o) are not allowed energetically. The isospin 
violating decay to Dgir together with electromagnetic decay to Dg'^ are then responsible for the 
narrow width observed. A similar situation may exist for Bg decays and we investigate this 
here using our lattice mass determinations of the excited states. This will enable us to predict 
whether narrow excited Bg mesons should be found. 

As well as exploring this issue of great interest to experiment, we determine the excited state 
spectrum of static-light mesons as fully as possible. This will help the construction of phe- 
nomenological models and will shed light on questions such as, whether there is an inversion 
of the level ordering with lighter than at larger I or for radial excitations as has been 
predicted [23l [Ml EH [Ml El] . 

Since we measure the spectrum for a range of values of the bare quark mass parameter /iq for the 
light quark, we could also compare with chiral effective Lagrangians appropriate to HQET. This 
comparison would be most appropriate applied to heavy-light decay constants in the continuum 
limit (see ref [28]). Since that study awaits more precise renormalization constants, we do not 
discuss it further here. 

Since we have discussed the basic methods in a previous paper [14j . in this paper we present 
only briefly the details of our computation of static-light meson mass differences. We give a full 
discussion of our extrapolation to the continuum and to physical light quark masses. We also 
discuss the interpolation to the physical b quark mass. 

2 Lattice details 

We use Nf = 2 flavour gauge configurations produced by the European Twisted Mass Collab- 
oration (ETMC). The gauge action is tree- level Symanzik improved [29] . while the fermionic 
action is Wilson twisted mass at maximal twist (cf. e.g. [30] and references therein). As argued 
in [13] this ensures automatic 0{a) improvement for static-light spectral quantities, e.g. mass 
differences of static-light mesons, the quantities we are focusing on in this work. 

We use three different values of the lattice spacing a ~ 0.051 fm , 0.064 fm , 0.080 fm and various 
values of the pion mass in the range 280 MeV < mps < 640 MeV. All lattice volumes are big 
enough to fulfill mpgL > 3.2. The ensembles we are considering are listed in Table [H Details 
regarding the generation of gauge configurations and analysis procedures for standard quantities 
(e.g. lattice spacing, pion mass) can be found in [STj [32]. 

In Table [T] we also list the number of gauges, on which we have computed static-light correlation 
functions, and the number and type of inversions performed to estimate light quark propagators 
stochastically. Note that in contrast to our previous work [131 [T4] we treat Bg mesons in a 
partially quenched approach, where the mass of the valence quark is approximately the mass of 
the physical s quark, taken from the study of strange mesons using the same configurations [HI 

US], 

• /3 = 3.90 /^q,valcncc = /^q.s = 0.022, 
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# and type 




L'-^ xT 




a in fm 


mps in MeV 


# of gauges 


of inversions 


3.90 


243 X 48 


0.0040 


0.0801(14) 


336(6) 


1420/580 


(spin) / (4rand) 






0.0064 




417(7) 


1480/- 


(spin)/- 






0.0085 




478(8) 


1360/480 


(spin) / (4rand) 






0.0100 




517(9) 


460/480 


(6rand) / (4rand) 






0.0150 




637(11) 


1000/- 


(Irand)/- 


4.05 


32^ X 64 


0.0030 


0.0638(10) 


321(5) 


240/240 


(4rand) / (4rand) 






0.0060 




443(7) 


500/500 


(4rand) /(4rand) 


4.20 


48^ X 96 


0.0020 


0.0514(8) 


284(5) 


420/420 


(spin) / (4rand) 



Table 1: ensembles (a and mp^ have been taken from [32]; ^ of gauges considered for B/Bg 
mesons; ^ and type of inversions for B/Bg mesons: (spin) four spin diluted timeslice sources 
on the same randomly chosen timeslice; (Irand) a single timeslice source on a randomly chosen 
timeslice; (4rand) four timeslice sources on four randomly chosen timeslices; (6rand) six timeslice 
sources on six randomly chosen timeslices). 



• /3 = 4.05 Aiq,valence = /^q,s = 0.017, 

• /3 = 4.20 /iq,valcncc = /^q,s = 0.015, 

while the sea is considerably lighter (cf. the listed //q values in Table H]) . 

3 Static-light mass differences 

The determination of static-light mass differences is essentially identical to what we have done 
in [laill]. 

For each of our ensembles characterised by the gauge coupling /? and the twisted light quark 
mass Hq (cf. Table [I]) and each of the lattice angular momentum representations Ai, E and A2 
we compute 6x6 static-light correlation matrices. The corresponding meson creation operators 
differ in their (twisted mass) parity, in their 7 matrix structure and in their spatial size. They 
are precisely the same we have been using before and are explained in detail in [H], section 3, 
Table 3. 

From these correlation matrices we compute effective mass plateaux using variational meth- 
ods [331 [M] (cf. |35j for exemplary plots showing the quality of our plateaus). We extract 
mass differences by fitting constants to these plateaus at sufficiently large temporal separations 
Tmin . . . Tmax- We determine Tmm and Tmax by requiring that the reduced is 0{1). T^nn 
values are listed in Table [21 while ?max — 

11 for /5 = 3.90 and /3 = 4.05 and T^ax = 17 for 
P = 4.20 in most cases (for some of the excited states smaller values had to be chosen, because 
the signal was lost in statistical noise). Note, however, that the choice of T^max 

is essentially 

irrelevant for the resulting mass (on the "Tmax side" of the effective mass plateau statistical 
errors are rather large and, therefore, data points only have a very weak effect on the fit). Since 
we are only interested in mass differences AM(j^) = M{j^) — M{S), the jackknife analysis has 
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been applied directly to the mass difference and not to the individual masses. The samples for 
M{j^) and M{S) entering for such a mass difference have been obtained with the same value 





B mesons 


Bs mesons 






P 


P+ 


D± 




F± 










D+ 


F± 


5* 


3.90 


0.0040 


6 


6 


5 


4 


4 


4 


6 


6 


5 


5 


4 


4 




0.0064 


6 


6 


5 


4 


4 


4 
















0.0085 


6 


6 


5 


4 


4 


4 


5 


5 


4 


4 


4 


4 




0.0100 


6 


6 


5 


4 


4 


4 


5 


5 


4 


4 


4 


4 




0.0150 


6 


6 


5 


4 


4 


4 














4.05 


0.0030 


7 


6 


6 


5 


5 


6 


7 


7 


6 


5 


5 


7 




0.0060 


7 


6 


6 


5 


5 


5 


7 


7 


6 


5 


5 


7 


4.20 


0.0020 


10 


10 


8 


7 


7 


9 


11 


9 


9 


8 


8 


11 



Table 2: Tmin for fitting constants to effective mass plateaus. 



The resulting mass differences AM(j^)a (in lattice units), where 

G {P- , P+ , D± , , F± , S*}, together with the pion masses mpsa (in lattice units; cf. 
Table [Hand [32]) and the lattice spacings a (in physical units; cf. Table [1]) serve as input for the 
extrapolation procedure to physical u/d quark masses described in the next section. 

We checked the stability of our results by varying Tj^in by il as well as by fitting superpositions 
of exponentials to the elements of the correlation matrices (as done in [1^) instead of solving 
generalised eigenvalue problems. We found consistency within statistical errors. 

4 Continuum limit and extrapolation to physical u/d quark 
masses 

4.1 Numerical results 

The mass differences AM(j^) obtained for all our ensembles are plotted against (mps)^ in 
Figure [U (unitary, i.e. mesons") and Figure [2] (partially quenched, i.e. "-B^ mesons"). Note 
that, although we use three different values of the lattice spacing, points corresponding to the 
same mass difference fall on a single curve. This is reassuring, since we use Wilson twisted mass 
lattice QCD at maximal twist, where static-light mass differences are 0(a) improved [13] • In 
Table [3] and Table [4] we collect the values of the mass differences in McVlj for all simulation 
points for B and Bg mesons respectively. 

For the extrapolation to physical light quark masses, we could use an effective field theory 
approach (Chiral HQET for instance) as used to study the decay constants [28] of the ground 
state. This approach has not been developed to discuss mass differences between excited states 

^The scale has been set by the pion decay constant as explained in detail in [32| . 
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(M{P_) - M{S))physi,a, = (406.03 ±18.74) MeV (x'^/d.o.f. = 0.95) 
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Figure 1: static- light mass differences linearly extrapolated to the physical u/d quark mass 
(unitary, i.e. B mesons). 
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(M{P_) - M{S))p^y,i,3| = (412.51 ± 1 1 .75) MeV (x'^/d.o.f. = 0.64) 



(M(PJ - M(S))physi,ai = (504.43 ± 1 1 .57) MeV (x^/d.o.f. = 0.64) 
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Figure 2: static-light mass differences linearly extrapolated to the physical u/d quark mass 
(partially quenched, i.e. Bg mesons). 
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/3 




AA/(P_) 


AA/(P+) 


AM{D±) 


AM{D+) 


AM{F±) 


AM{S*) 


3.90 


0.0040 


415(17) 


494(20) 


855(30) 


879(25) 


1155(35) 


749(22) 




0.0064 


449(17) 


499(20) 


879(26) 


924(24) 


1253(33) 


740(21) 




0.0085 


471(17) 


506(19) 


878(25) 


928(24) 


1223(40) 


766(20) 




0.0100 


474(22) 


481(21) 


881(34) 


889(32) 


1225(40) 


755(23) 




0.0150 


513(29) 


465(21) 


829(50) 


889(48) 


1192(45) 


794(24) 


4.05 


0.0030 


465(26) 


495(24) 


887(39) 


952(49) 


1148(60) 


821(44) 




0.0060 


498(22) 


551(23) 


851(44) 


1000(41) 


1273(53) 


794(20) 


4.2 


0.0020 


399(31) 


498(35) 


851(45) 


990(53) 


1184(58) 


845(51) 



Table 3: static-light mass differences in MeV (unitary, i.e. B mesons) for all simulation points; 
details on the analysis procedure of the correlation functions are given in section [3j 







AM(P_) 


AM(P+) 


AM{D±) 


AM(P>+) 


AA/(P±) 


AM{S*) 


3.90 


0.0040 
0.0085 
0.0100 


438(13) 
466(14) 
471(15) 


499(14) 
495(14) 
497(13) 


805(30) 
888(23) 
882(20) 


902(35) 
880(24) 
855(28) 


1193(37) 
1171(41) 
1219(40) 


729(26) 
730(21) 
726(22) 


4.05 


0.0030 
0.0060 


444(13) 
422(14) 


500(13) 
491(13) 


810(26) 
842(23) 


934(24) 
918(22) 


1167(36) 
1223(32) 


734(29) 
735(31) 


4.2 


0.0020 


417(13) 


509(13) 


811(29) 


930(34) 


1226(52) 


790(41) 



Table 4: static-light mass differences in MeV (partially quenched, i.e. Bg mesons) for all simula- 
tion points; details on the analysis procedure of the correlation functions are given in section [3l 



and the ground state (e.g. M(P_) — M{S)), so is not appropriate here. Instead we use the 
simplest assumption which is supported by our results: a linear dependence. 

Because our ground state mass values enter into all of the mass differences we study, we simul- 
taneously fit to all the meson mass differences we have computed. We find that fits which are 
independent of the lattice spacing and which are linear in the light quark mass (represented by 
the mass squared of the light-light pseudoscalar meson) are acceptable, i.e. yield x^/dof<l. 

For the Bs mesons, our results depend on the strange quark mass we choose. We have taken these 
values from studies of strange-light mesons [15\ [IB] as discussed above. The possible systematic 
error arising from an incorrect value for the strange quark mass is very small: because the mass 
differences we measure turn out to be very weakly dependent on that mass. This will be seen 
when we compare our results for the B and Bg mesons extrapolated to physical light quark 
masses. 

The details of our fitting procedure are collected in appendix [XI 

As already mentioned both fits (one for B mesons, the other for Bg mesons) are of good quality in 
a sense that x'^/dof < 1. This shows that at the present level of statistical accuracy the continuum 
limit has already been reached at our largest value of the lattice spacing a ~ 0.080 fm. Moreover, 
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these fits enable us to extrapolate to physical u/d quark masses. 

Extrapolations of static-light mass differences to physical u/d quark masses are listed in Table [5] 
in MeV both for B mesons and for Bg mesons. Note that both fits give x^/d.o.f. ~ 1, i.e. are 
consistent with our assumption that static-light meson mass differences as functions of (mps)^ 
can be parameterised by straight lines. 





P_ 


P+ 


D± 


D+ 


F± 


5* 


xVd.o.f. 


B mesons 


406(19) 


516(18) 


870(27) 


930(28) 


1196(30) 


755(16) 


0.95 


Bg mesons 


413(12) 


504(12) 


770(26) 


960(24) 


1179(37) 


751(26) 


0.64 



Table 5: M{y) — M{S) in MeV extrapolated to physical light quark masses. 



To check the stability of these fits, we have varied T^\^ by ±1. Within statistical errors mass 
differences obtained with Tmin — 1, with Tmin and with Tmin + 1 are in agreement. 

The extrapolations are shown in Figure [1] and Figure [2j The red dots represent the maximum 
likelihood estimates of z = ((mps)^ , AM(j^)) obtained during the fitting procedure. In ad- 
dition to x-y-error bars we also plot covariance ellipses, which reflect the correlations between 
(?7T,ps)^ and AM(j^) induced by the lattice spacing a, that is they are generated from the 
inverses of the corresponding 2x2 submatrices of the covariance matrix C. 

4.2 Contamination of static-light meson masses by multi particle states 

The radially and orbitally excited static-light mesons P_, P4., Z?+, F„, -F+ and S* can 
decay into multi particle states S" -|- n x vr with relative angular momentum such that quantum 
numbers are identical. In particular the P_ static-light meson is not protected by angular 
momentum, i.e. it can decay via an S wave into 5 -|- vr, whose wave function is not suppressed 
at the origin. In the following we argue that the effect oi S -\- tt states on our P„ mass is small 
compared to its statistical error. To this end we resort to a model presented and to numerical 
results obtained in [36 | [37 l [38]. 

We consider the P_ static-light meson at /3 = 3.90 and our lightest u/d quark mass at this (3 
value (^q = 0.0040). In that ensemble the masses of the P_ state and of the S + t: state are quite 
similar: m[P^)a ~ 0.57 and {m[S) + m(7r))a ~ 0.53 (we consider the case, where the pion has 
zero momentum). Therefore, we expect mixing of P_ and S + 1: with respect to the eigenstates 
of the Hamiltonian H, mixing which will be different in different spatial volumes. Consequently, 
we do not focus on the eigenvalues of these states, but rather on m(P_) = {P-\H\P-) {\P-) 
is a state with = (1/2)+ created by single particle operators, e.g. operators of type Qu or 
Qd, which we have used in the construction of trial states). At very large temporal separation 
the correlators we are studying will inevitably yield the eigenvalues of the Hamiltonian. At 
intermediate temporal separations, however, one can expect to read off m(P_) as we will explain 
in the following. 

In [38] the effective coupling strength of the decay P_ — )■ 5 -|- vr has been estimated by a lattice 
computation: F/Zc ~ 0.46. Moreover, some evidence has been obtained that this quantity is 
fairly independent of the light quark mass. Using this result one can determine the mixing 
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element xa of the energy matrix via eqn. (5) in |38] for our situation {L/a = 24, iJij^a ~ 0.14): 

/ 2TT(r/k) 

\3{L/ay{m^a)J 

As detailed in |36t [571 ES] for large temporal separations the P_ correlator is of the form 

Cp_{t/a) cx e-^'"^- cosh((xa)(Va)), (2) 
while the corresponding effective mass is 

"T'effcctivc,p_(Va)a = - -^^jjj-^ In (^Cp_{t/a)^ = 



d 



(mp_a)(t/a) — In ^cosh((xa)(t/a))^ j = mp_a — tanh((xa)(t/a))xa. (3) 



d{t/a) 

At t/a = 12 (the maximum temporal separation we have considered) the estimated systematic 
error of mp_ coming from mixing with S" + tt is tanh((xa)(t/a))xa ~ 0.0063, i.e. roughly a 1% 
effect. This correction is significantly smaller than the statistical error of meffective,P_o in that t 
region. 

For the other temporal separations and/or ensembles we obtain similar estimates. We, therefore, 
expect that at the present level of statistical accuracy the effect of multi particle states on our 
static-light meson masses, in particular on P_, is negligible. 

Our conclusions are in agreement with those obtained in [9], where a study of the static-light 
meson spectrum with similar techniques has been performed using two different lattice volumes. 
No volume dependence of the eigenvectors of static-light meson states has been observed, which 
is a sign that contributions of multi particle states are strongly suppressed. 



5 Extrapolation to the physical b quark mass 

To make contact with experimentally available results on the spectrum of B mesons, we need to 
correct for the non-infinite mass of the b quark. In Heavy Quark Effective Theory, the leading 
correction will be of order 1/niH, where mn is the heavy quark mass. It is possible, in principle, 
to evaluate the coefficients of this correction from first principles on a lattice [191 120] . This we 
intend to explore in the future, but here we use a more direct method to establish the size of 
this small correction between static quarks and b quarks of realistic mass. These l/mn terms 
will break the degeneracy of mesonic states found in the static limit. 

We evaluate for physical b quarks by interpolating between static heavy quarks and the charm 
quark, where experimental data is available. As a measure of the heavy quark mass, we take the 
mass of the ground state heavy-light meson (D or B). This measure is equivalent to another 
(such as using quark masses in some scheme) to the order l/mn we are using. One test of 
this interpolation can be made. The hyperfine splitting between D* and D of 141 MeV when 
interpolated from the static limit (namely zero) gives for B* and B a splitting reduced by 
m{D)/m{B) = 0.35 to 49 MeV which agrees with the observed splitting [22] of 46 MeV to within 
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For the fine splitting, the kinetic term (rather than the chromo-magnetic) is relevant and the 
experimental results for the spectrum are rather incomplete - indeed this current study is to 
establish the spectrum from a theoretical input. Lattice studies do confirm |19|l2Uj that a l/mn 
behavior is dominant down to masses near the charm quark mass. 

We interpolate our lattice results for static-light mass differences of P and S wave states to the 
physical b quark mass at m{D)/m[B) = 0.35 linearly in m{D)/mH, making use of experimental 
data on D and Ds mesons as input [22j . For details regarding this method of extrapolation cf . 
|14| . Results are listed and compared to experimental results in Table O The corresponding 
extrapolations are shown in Figure [3j 

For D mesons the assignment of the two = l"*" states to i?^ and Bi is easy, because their 
widths differ by more than an order of magnitude (we associate the narrow state with Bi [one 
of the two degenerate = (3/2)"^ states in the static limit, which can only decay to S' + vr via a 
D wave and is, therefore, protected by angular momentum]; the wide state with B* [one of the 
two degenerate = (1/2)"'' states in the static limit, which can readily decay to 5 + 7r via an S 
wave]). In contrast to that the situation is less clear for Ds mesons, where both = l"*" states 
have similar (narrow) widths. Therefore, we show both possibilities in Tableland in Figure [3l 





M - M{B) in MeV 




M - M{Bs) in MeV 


state 


lattice 


experiment 


state 


lattice 


experiment 


B* 


443(21) 




BsO 


391(8) 




BI 


460(22) 




B*si 


440(8)/467(8) 




Bi 


530(12) 


444(2) 


Bsi 


526(8)/499(8) 


463(1) 


B*2 


543(12) 


464(5) 


Bs2 


539(8) 


473(1) 


B* 




418(8) 


Bh 




487(15) 



Table 6: lattice and experimental results for P wave B and Bg states {Bj and B*j denote 
rather vague experimental signals, which can be interpreted as stemming from several broad 
and narrow resonances possibly including the j = 1/2 P wave states Bq, B\, B*q and B*^] 
the two lattice values listed for B^^ and Bsi correspond to the two possibilities of assigning 
experimental = D results [cf. text for more details]). 

Compared to our previous study [131 at a single lattice spacing, the above results are similar 
for the B (unitary) case. For Bs mesons we now employ a partially quenched s quark which 
allows a more realistic treatment of the light quark sea. So our new results supersede those 
obtained previously for Bs- Indeed we find a significant dependence on the sea quark mass (cf. 
Figure [2]), which is now the physical u/d quark mass, while it previously corresponded to the 
significantly heavier s quark mass. 

In our lattice study we have extracted the continuum limit and have extrapolated to physical 
light quarks using a linear dependence. We have then interpolated to the physical h quark 
assuming that a l/mn behavior is valid down to the charm quark mass. These assumptions 
induce systematic errors and, in principle, they can be quantified by further lattice studies. 

The assumption of a linear extrapolation to physical light quarks is sensitive to possible admix- 
tures of two body states which become more important at lighter quark masses as thresholds 
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a) prediction tor excited B states 
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Figure 3: Static-light mass differences linearly extrapolated to the physical b quark mass, a) Uni- 
tary, i.e. B mesons, b), c) Partially quenched, i.e. Bg mesons. 
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for decay open. We have explored this possibility and found no evidence of such effects, so it 
is difficult to estimate the magnitude of a possible systematic error from this. If there was a 
significant difference between the light quark behavior for the ground state and an excited state, 
this would introduce an error on our extrapolation to the physical value which could be as large 
as lOMeV. 

The test of the l/mn assumption for the chromo-magnetic term, discussed above, was found to 
be valid within 6%. This suggests that an estimate of the systematic errors for the B and Bg 
meson mass splittings coming from 1/m effects should also be at least of order 6%. Since the 
1/m correction to the P wave states is of order 100 MeV, this implies a systematic error of order 
6MeV. 

One further possible source of systematic error is from our neglect of the strange contribution 
to the sea. This will be addressed in a future study making use of the Nf = 2 + 1 + 1 sea which 
includes dynamical s quarks from ETMC \39\ HQ] . 

Overall, it seems prudent to assign systematic errors on our mass differences (for P_ and 
relative to S) of order 20 MeV from these effects, even though we have little evidence for such 
effects. 

The experimental determination of the spectrum of excited B and Bg mesons is quite limited |22j . 
Assuming that the relatively narrow states seen correspond to our P+ state (since a = 
state must have that assignment), the mass difference we see of over 500 MeV does not agree 
closely with the experimental results of around 450 MeV. We do get a mass difference of around 
450 MeV from our P_ states, although such states cannot have = 2+. 

In view of this discrepancy with experimental results, it is also interesting to compare with 
independent existing lattice computations, in particular with the rather recent study reported in 
|12j . There the light quark extrapolation is only performed in the valence quark mass (from which 
static-light mass differences essentially seem to be independent, as can be seen by comparing 
our B and Bg results and also from corresponding plots and numbers presented in [12j). while 
the sea quark mass is kept fixed. More generally, a comparison of the dependence of static-light 
mass differences on the sea quark mass, which we have computed down to mps ~ 280 MeV, 
with existing lattice studies is not possible: there the number of investigated sea quark masses 
is rather small and they are quite heavy, around the mass of the s quark. What one can do, 
however, is to compare meson mass differences for a given value of the sea quark mass. Before 
comparing results (in physical units) with those quoted in |12j it should be noted that in jl2j the 
scale is set by identifying tq with 0.49 fm, while our result for this quantity is ro = 0.42 fm [32) . 
Therefore, to perform a meaningful comparison, one should express all quantities in units of ro or 
equivalently scale all masses in physical units listed in [12] by a factor of around 0.49/0.42 ~ 1.14. 
For the lightest sea quark mass considered in |12] corresponding to mpg ~ 461 MeV it is most 
appropriate to compare with our results at /3 = 3.90, = 0.0100 (mps 517 MeV). For the P 
wave mass differences one finds 

{m{P_) - m(5'))ETMC 

(m(P_)-m(5))[T2] 
(m(P+) - m(5'))ETMC 
(m(P+) -m{S))^ 

ratios, which are within statistical errors fully consistent with the expected factor 1.14. 



454(19)(9)MeV -^-^^1^-^^ l^J 

481(27) MeV 

^ ^ - 1.08(11), (5) 



446(17) (9) MeV 
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It is interesting to note that the ratios between our lattice results and the experimental values (see 
Table [6]) are on the same ballpark of the ratio between two values of ro used above, i.e. ~ 1.14. 
While there is no reason to doubt the precise determination of the lattice spacing performed in 
[32], it would be interesting, although beyond the scope of this paper, to investigate, whether 
simulations at lighter quark masses and /or with Nf = 2 + 1 + 1 dynamical flavours will improve 
the agreement with experimental results. 

One interesting issue is whether the Bg states are stable to the strong decay to BK. This decay 
has a threshold at 408 MeV above the ground state Bg meson. Our states (the upper two in 
Table [6|) do indeed have masses which are close to (or below) this threshold. That would imply 
that these two states {B*q and B*^) should have a very small decay width. This is consistent 
with the experimental observation that only two candidate P wave B* states have been seen so 
far: corresponding to the heavier P^ states. All the other states Bg we study, including the S*, 
lie higher than this BK threshold and so would have a strong decay open. 

Moreover, our findings clearly indicate that there is no inversion of level ordering for P wave 
states, neither for B mesons nor for Bg mesons. Bq and B^ {B*q and -B*^) are considerably 
lighter than Bi and B2 (Bgi and -B52) ^™ be read off from Tableland Figure [3l This is in 
contrast to predictions obtained from certain phenomenological models [231 [24} [25| [26| [27] and, 
therefore, might provide valuable input for future model building. 

6 Conclusions 

We have determined the continuum limit for static-light mesons on a lattice using Nf = 2 
flavours of light quarks. The removal of 0{a) effects by using maximally- twisted mass fermions 
for meson mass differences in the static limit is confirmed. 

We have investigated the light sea quark mass dependence of B and Bg mesons down to 
mps ~ 280 MeV, which is significantly lighter than what has been achieved in previous studies 
of static-light mesons. We find that our results are compatible with a linear extrapolation in the 
light quark mass to its physical value. We see no sign of any mixing with two body effects and 
this is consistent with our estimate that such effects should be too small to see on our lattices. 

We have determined masses for a wide variety of excited states in the continuum limit and this 
will be a valuable resource for model builders. 

We have employed the assumption of a l/rriH dependence on the heavy quark mass together 
with experimental results for charm-light mesons to allow us to estimate the spectrum that one 
would obtain for physical b quarks. 

Our results imply that there will be a = 0^ and = 1+ Bg meson which has a narrow 
width since its strong decay to BK is suppressed (or zero) due to phase space effects. 

Future directions include (i) determination of Jb and Jb^ (for a preliminary result cf. [28]): 
(ii) a similar investigation regarding static-light baryons; (iii) extending these computations to 
Nf = 2 + 1 + 1 flavour ETMC gauge configurations [39l[4n] . 
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A Details of the fitting procedure 



Data points (mps)^ and AM{j'^), G {P_ , P+ , D± , , F± , S*} corresponding to the same 
/3 are correlated via the lattice spacing a. We take that into account via a covariance matrix, 
which we estimate by resampling nipsa, /S.M{j^)a and a (100,000 samples). Consequently, we 
do not fit straight lines to the Data points ((mps)^ , AM(j^)) individually for every static-light 
state , but perform a single correlated fit of six straight lines to the six mass differences of 
interest. During the fitting we take statistical errors both along the horizontal axis (errors in 
(mps)^) and along the vertical axis (errors in AM(j^)) into account. 

The method of performing the two-dimensional fits is based on what has been used in |41| . 

To be able to express the corresponding equations in a compact way, we introduce the following 
notation: 

. z = (x, y(l), y(2), ...). 

• X = (((w-ps)^)^"*^^ , ((w-ps)^)^^^ 1 • • •) (the upper index refers to both the lattice spacing 
and to the light quark mass). 

• y(j) = ((AM)(i)(j), (AM)(2)(j), ...) (the upper index (■■■) refers to both the lattice 
spacing and to the light quark mass, the index (j) refers to j^). 

• C denotes the estimated covariance matrix for z (a 56 x 56 matrix for B mesons, a 42 x 42 
matrix for Bg mesons). 

• The linear fits y{j) = a{j)x + b{j) are parameterised by a(j) and b{j) (the quantities, 
which will finally allow the extrapolation to physical u/d quark masses). 

The basic idea of the method is a maximum likelihood determination of the "true values" 
z = (x , y(l) , y(2) , . . .). This amounts to minimizing 

1 / \T 



— z — z 

2 



C-^(z-z) -^A„(i)(a(i)x„ + 6(i)-y„(i)) (6) 



with respect to z, a(j), b{j) and A(j) under the constraints Unij) = CL{j)xn + b{j). For z we 
use the same resampling procedure as for estimating the covariance matrix (this is necessary, 
because za = (((mps)^)^) / ((mps)(")a)V(a)' and za = ((AM)(")(j)) / ((AM)(")(j»/(a)). 

The constraint minimization is equivalent to solving a system of non-linear equations, which we 
do by means of the scaled- hybrid algorithm of the GSL library |42j . It needs initial parameters, 
which should preferably be close to the global extremum. Such initial parameters can be obtained 
by individual standard one-dimensional straight line fits: 

• A„(i) = 0, 

• a{j) and b[i) minimizing 
„ {a{j)xn + b{j) - yn{j)) 

r 
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• X = X and y{j) = y{j). 



To judge the quality of the resulting fit, we define a "reduced 



via 





(8) 



d.o.f. 



where d.o.f. is the number of entries of z minus the number of a{j) and b{j), i.e. d.o.f. = 44 for 
B mesons and d.o.f. = 32 for Bg mesons respectively. 

The resulting straight lines allow an extrapolation to physical u/d quark masses (corresponding 
to mps = 135 MeV). The corresponding statistical errors are obtained by repeating this fitting 
and extrapolation procedure 100 times with randomly sampled sets za (we randomly sample 
the input data and compute za = ((mps)'^")a)^/a^ and za = {{AM)^"'\j)a) / a) and taking the 
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